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Materials with zero/near zero thermal expansion coefficients are technologically important for applications 
in thermal management and engineering. To date, this class of materials can only be produced by chemical 
routes, either by changing chemical compositions or by composting materials with positive and negative 
thermal expansion. Here, we report for the first time a physical route to achieve near zero thermal expansion 
through application of pressure. In the stability field of tetragonal PbTi0 3 we observed pressure-induced 
reversals between thermal contraction and expansion between ambient pressure and 0.9 GPa. This hybrid 
behavior leads to a mathematically infinite number of crossover points in the pressure-volume-temperature 
space and near-zero thermal expansion coefficients comparable to or even smaller than those attained by 
chemical routes. The observed pressures for this unusual phenomenon are within a small range of 0.1- 
0.9 GPa, potentially feasible for designing stress-engineered materials, such as thin films and nano-crystals, 
for thermal management applications. 



Most materials expand on heating and contract on cooling. Although uncommon, some materials exhibit 
opposite behavior and possess negative coefficients of thermal expansion (CTE). There is yet a third 
family of materials that neither expand nor contract as temperature changes, and they display zero 
thermal expansion (ZTE) within a certain temperature range 1 " 6 . ZTE is an intriguing and useful physical attribute, 
because materials with ZTE do not undergo thermal shock or fatigue during rapid or repeated heating and cooling 
cycles, which makes them excellent candidates for applications in optics, electronics, and heat-engine compo- 
nents. ZTE is typically achieved by composting materials with positive and negative thermal expansion (NTE), 
but this route is often hindered by poor thermal stability of materials with NTE 7 . The mismatch in thermal 
expansion also results in internal strains in composite materials at elevated temperatures. Chemical substitution, 
operating as "chemical stress", is another common method to modify structural and physical properties of 
materials 8 ' 9 . In AB0 3 perovskite, for example, chemical substitutions over both the A and B sites are effective 
pathways to tune materials with NTE or even ZTE over a large temperature range, as demonstrated by the studies 
in PbTi0 3 -based pervskites 10 " 15 . ZTE phenomenon is also fundamentally interesting because it is associated with 
rich and complex physical mechanisms 12 ' 4 " 6 , such as Invar effect, valence-state transition, strain relaxation of the 
constituent polyhedra in crystal structures, and elastic anisotropy. 

Lead titanate, PbTi0 3 , is one of the most widely used ferroelectric and piezoelectric materials in industry 16 " 18 . 
At ambient conditions, PbTi0 3 possesses a tetragonal structure (P4mm) with an axial ratio of c/a= 1.06 19 . The 
structural distortion from the ideal, cubic structure (space group Pm3m) is generally attributed to the covalent 
nature of B-O bonds in AB0 3 -type ferroelectric perovskite. In tetragonal PbTi0 3 , the 6s 2 lone pairs effect of Pb 2+ 
and the Pb-O hybridization are not only responsible for the extra degree of structural distortion 20 ' 21 , high 
spontaneous polarization, and relatively high Curie temperature (763 K) 22 , they are also important mechanisms 
underlying the observed NTE 10 " 14 ' 23 . The mean CTE of this compound is — 1.99 X 10~ 5 K" 1 over the temperature 
range of 300-763 K 15 , and the corresponding thermal contraction is closely related to the ferroelectric soft modes 
with increasing temperature 15 . Above 763 K, the tetragonal PbTi0 3 transforms to a paraelectric cubic phase that 
exhibits the regular positive CTE 17,24 " 28 . 
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Figure 1 | A representative neutron diffraction pattern of PbTi0 3 
collected at 480 K and 0.71 GPa. The observed and refined profiles are 
denoted, respectively, by open circles and red lines. NaCl is used for 
pressure calibration and C (graphite) is from the furnace. 

As an important thermodynamic variable, pressure can effectively 
change the bonding characteristics (bond lengths and angles) and 
electron clouds overlapping of the constituent atoms in a compound, 
which would ultimately alter the material's structural and physical 
properties. In PbTi0 3 , both tetragonality and ferroelectricity are 
progressively suppressed with increasing pressure at 300 K, leading 
to a paraelectric cubic phase at pressures above 12 GPa 29 . Compared 
with chemical substitution, the application of pressure is a straight- 
forward route to tune material properties without introducing 
impurities and/or vacancies into the parent structure. Such defects 
are known to affect many properties such as structural stability, 
chemical diffusion, electric and ionic conductivities, and elasticity. 
In the present work, we applied external pressure to manipulate 
CTE of tetragonal PbTi0 3 , including the achievement of near zero 
thermal expansion. This approach was motivated by our earlier find- 
ings 30 that at pressures above 1 GPa PbTi0 3 expands at all tempera- 
tures within the stability field of tetragonal phase. This behavior 
implies a crossover in the volumetric thermal expansion from nega- 
tive to positive between ambient pressure and 1 GPa, which would 
presumably lead to mean ZTE over certain temperature ranges. Such 
pressure-tuned ZTE, if it can experimentally be verified, is not only 
scientifically important but also technologically relevant because 
similar levels of stress can readily be realized in thin-film and nano- 
crystalline PbTi0 3 produced by stress engineering 31 " 33 . 

Results 

Figure 1 shows a representative neutron diffraction pattern of 
PbTi0 3 analyzed by the Rietveld method. The refined structural 
parameters, such as lattice parameters, atomic positions, goodness 
of fit and R factors, are listed in Table 1. Figure 2a shows the pressure 
dependence of unit-cell volumes of tetragonal PbTi0 3 at different 
temperatures. In striking contrast to the compression behavior in 
most materials, we observe a number of crossovers between the P- 
V isotherms within the 0-1 GPa range (Fig. 2a), implying reversals in 
the sign of CTE of PbTi0 3 . At 0.1 GPa, there is a crossover point 
between the 300 K and 480 K volume isotherms, indicating that the 
mean volumetric CTE, a v = (V 2 - V 1 )/V 1 (T 2 - Tj), is zero over this 
temperature range. Similar to the behavior at ambient pressure, tet- 
ragonal PbTi0 3 shrinks at temperatures above 480 K with a negative 
CTE. The second crossover point is located at —0.2 GPa between the 
300 K and 580 K isotherms. PbTi0 3 , however, expands with increas- 
ing temperature between 300 and 480 K and continues to contract at 



higher temperatures. The third crossover point is found at 
—0.55 GPa between the 300 and 730 K isotherms. At this pressure, 
the structure expands between 300 and 580 K with essentially zero 
thermal expansion between 480 and 580 K isotherms; the ot v is zero 
over the 300-730 K range due to thermal contraction between 580 
and 730 K. Evidently, at a given pressure, tetragonal PbTi0 3 exhibits 
a hybrid behavior of contraction and expansion. Such unusual beha- 
vior eventually leads to three additional crossovers at higher pres- 
sures, 0.65 GPa between 480 K and 580 K isotherms, 0.85 GPa 
between 480 K and 730 K isotherms, and 0.9 GPa between 580 K 
and 730 K isotherms. When pressure is above 0.9 GPa, the unit-cell 
volumes expand at all experimental temperatures within the stability 
field of tetragonal PbTi0 3 , which is consistent with our earlier obser- 
vations using synchrotron XRD 30 . 

As a result of pressure-induced reversal between thermal contrac- 
tion and expansion, there are mathematically an infinite number of 
crossover points in the P-V-T space, giving the infinitesimal ATs 
between two volume isotherms, where AT = T 2 - 7\ with T 2 > 7\. 
Our findings indicate that the larger the mean temperature of T\ and 
T 2 , the higher the pressure is needed to produce the crossover 
between the corresponding isotherms (Fig. 2a). This behavior would 
in principle lead to mean ZTE coefficients over an infinite number of 
ATs. The exact temperature range for ZTE or near ZTE, however, 
cannot be determined based on the measured data because pressure 
varies with heating in complicated manners that depend on thermal 
pressure, deviatoric stress, and materials' strength. To determine the 
temperature ranges for ZTE, we interpolated the experimental data 
of Figure 2a and plot the lattice volumes as a function of temperature 
in Figure 2b at three constant pressures of 0.3, 0.8, and 1.2 GPa, 
arbitrarily chosen for clarity. Strictly speaking, there are no such 
temperature ranges in which PbTi0 3 exhibits truly zero thermal 
expansion (except at the "turnover" points of 447 and 556 K where 
CTE changes from positive to negative). In other words, CTEs would 
deviate from ZTE when the temperature is away from the "turnover" 
points. However, pressure plays a significant role in changing ther- 
mal behavior from contraction at ambient pressure to expansion 
above 0.9 GPa. In addition, pressure tends to elevate the temperature 
of turnover point, which leads to increased temperature ranges for 
the expansion (AT Ei ) and reduced ranges for the contraction (AT Ci ) 
as pressure increases. Consequently, within the 150 K range of the 
turnover points, the lattice volumes of PbTi0 3 are nearly independ- 
ent of temperature. Based on such interpolated data, we calculated oc v 
values at pressures of 0.3 GPa and 0.8 GPa, which are listed in table 2 
along with those previously obtained from chemical substitutions for 
PbTi0 3 -based perovskite. Strictly speaking, neither chemical nor 
pressure-tuned route produces thermal expansion coefficients that 
are truly zero. Nevertheless, the resultant CTEs are substantially 
smaller than those for chemically pure PbTi0 3 at ambient pressure. 
At both pressures, the pressure-tuned "near-zero" CTEs are com- 
parable to or even smaller than those attained by chemical routes, 
indicating that pressure-tuning can be at least as effective as chemical 
routes for controlling materials' thermal behavior. 

Discussion 

In most materials, thermal expansion decreases with increasing pres- 
sure, which would in principle lead to a vanishing a at a finite pres- 
sure 34 . However, this general phenomenon does not allow the 
crossovers in pressure-volume isotherms. In fact, the extrapolation 
of experimental P-V-T data often requires some additional empirical 
constants to prevent such crossovers from occurring 35 . In the present 
work, we demonstrate for the first time that the P- V isotherms can 
crossover in materials with NTE at ambient pressure. The pressure- 
induced reversal from thermal contraction to expansion also repre- 
sents a new physical route for tuning near- zero CTEs. 

The pressure-tuned near-zero CTE in PbTi0 3 is technologically 
relevant to applications in thermal management and engineering. 
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Table 1 | The refined structural parameters of tetragonal PbTi03 and refinement agreement 

Tetragonal PbTi0 3 o o o Ti/O atom position (Z) Pb (0, 0, 0), Ti ( 1 /2, 

Pressure/Temperature a (A) c (A) V (A 3 ) 1 /2, Z) Ol (1 /2, 1 /2, Z) Oil (1/2, 0, Z) Refinement agreement indices 



1 atm/RT 


3.9154(2) 


4.1622(7) 


63.810(1) 


0.226 GPa/RT 


3.9214(6) 


4.1321(8) 


63.541(2) 


0.652 GPa/RT 


3.9189(7) 


4.1173(10) 


63.231(3) 


0.760 GPa/RT 


3.9185(8) 


4.1 141(1 1) 


63.171(3) 


0.858 GPa/RT 


3.9180(9) 


4.1093(11) 


63.083(3) 


1.21 1 GPa/RT 


3.9170(9) 


4.0968(11) 


62.857(4) 


1 atm/480 K 


3.9360(4) 


4.1 155(9) 


63.758(1) 


0.352 GPa/480 K 


3.9315(3) 


4.1 142(9) 


63.593(12) 


0.710GPa/480 K 


3.9319(3) 


4.0962(10) 


63.326(13) 


0.906 GPa/480 K 


3.9282(4) 


4.0965(1 1) 


63.211(14) 


1 .027 GPa/480 K 


3.9281(4) 


4.0926(12) 


63.149(15) 


1 .297 GPa/480 K 


3.9281(4) 


4.0792(13) 


62.943(17) 


1 atm/580 K 


3.9459(1) 


4.0905(3) 


63.690(4) 


0.434 GPa/580 K 


3.9398(3) 


4.0902(9) 


63.490(12) 


0.784 GPa/580 K 


3.9376(3) 


4.0824(10) 


63.297(13) 


0.940 GPa/580 K 


3.9390(4) 


4.0741(13) 


63.213(16) 


1 .080 GPa/580 K 


3.9366(4) 


4.0745(12) 


63.141(15) 


1.361 GPa/580 K 


3.9379(5) 


4.0613(14) 


62.979(18) 


1 atm/730 K 


3.9580(2) 


4.0534(37) 


63.501(2) 


0.438 GPa/730 K 


3.9557(8) 


4.0462(9) 


63.313(12) 


0.856 GPa/730 K 


3.9547(3) 


4.0441(10) 


63.246(13) 


0.978 GPa/730 K 


3.9556(4) 


4.0409(10) 


63.227(14) 


1.109 GPa/730 K 


3.9554(4) 


4.0405(1 1) 


63.214(14) 


1.541 GPa/730 K 


3.9591(6) 


4.0456(15) 


63.126(18) 



The required pressures for these phenomena are in a small range of 
0.1-0.9 GPa, which can readily be realized in thin films and nano- 
crystals manufactured by stress -engineering processes. It is well 
known that the residual stress in thin films can be manipulated by 
their thickness or mismatch with substrate materials. For PbTi0 3 
deposited on a Si(100) substrate 31 , the stresses in thin films of 
407 nm and 92 nm thickness are 0.96 (9) GPa and 1.77(12) GPa, 
respectively. They can be increased to the 1.3-2.6 GPa range when 
PbTi0 3 thin films of 400-50 nm thickness are deposited on a Pt- 
coated Si substrate 32 . It is well-known that in thin films the in-plane 
lattice parameters are constrained by the substrate materials. 
However, as a result of Pb-O and Ti-O bonding hybridizations, the 
a and c axis are coupled in the PbTi0 3 thin film, as demonstrated in 
ref [36]. Therefore, the stress can still affect the out-of-plane lattice 
parameter and manipulate the CTEs of PbTi0 3 thin films. Nano- 
synthesis can also introduce different levels of residual stress on 
the surface layer of a nanocrystal. The estimated stress in PbTi0 3 
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nanocrystals of 15 nm can be as high as 0.9 GPa based on the com- 
parison of cla ratios between this work and Ref. [33]. Thus, the 
present findings not only offer a physical route to tune near- zero 
CTEs in chemically-pure PbTi0 3 but may also potentially be useful 
for designing stress -engineered materials, such as thin films and 
nano-crystals, of desired thermal expansion coefficients, including 
ZTE. Last but not least, we expect that pressure-induced ZTE may 
occur in other NTE materials if they are structurally stable over the 
desired P-T space. 

The thermal expansion behavior and associated ferroelectric 
polarization in PbTi0 3 perovskite are commonly linked with its 
tetragonality and chemical bonding. To understand the correlation 
between these characteristics and the present findings in the pres- 
sure-temperature space unexplored by previous studies, we look into 
the details of cla ratios and hybridization between Pb and O atoms. 
Plotted in Figure 3a are the lattice parameters as a function of pres- 
sure for tetragonal PbTiQ 3 at 300 K, 480 K, 580 K, and 730 K. With 
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Figure 2 | (a) The pressure dependence of unit cell volumes of tetragonal 
PbTi0 3 at different temperatures. (2b), the lattice volumes of tetragonal 
PbTi0 3 as a function of temperature at constant pressures. The solid curves 
at 0.3 GPa, 0.8 GPa and 1.2 GPa are the data interpolated from Figure 2a. 
The plotted symbols are the experimental data points collected within 
±0.1 GPa of the indicated pressures. The pink shadow area denotes the 
region of thermal contraction, and the rest plain area corresponds to the 
region of thermal expansion. The AT Ei represents the temperature range 
for expansion and AT ci the temperature range for contraction. The 
subscripts of 1, 2, 3, and 4 in AT Ei and AT ci correspond to the pressures of 
0 GPa, 0.3 GPa, 0.8 GPa and 1.2 GPa, respectively. The turnover points 
from expansion to contraction are 447 K at 0.3 GPa and 556 K at 0.8 GPa. 



increasing pressure, both lattice parameters a and c decrease 
(Figure 3a). However, the rate of decrease in c is considerably larger 
than that in a, leading to smaller cla ratios at high pressures. 
Consequently, as the pressure increases, the tetragonal PbTi0 3 
approaches the cubic structure. Figure 3a also reveals the effect of 
temperature on the structural distortion; at a given pressure, lattice 
parameter a increases with increasing temperature whereas c 
decreases. In addition, the magnitudes of these variations decrease 
with increasing pressure. Figure 3b shows that the cla ratios are 
suppressed by both temperature and pressure, indicating that pres- 
sure can effectively modulate the tetragonality of PbTi0 3 . 

The thermal expansion behavior in tetragonal PbTi0 3 under 
atmospheric pressure is commonly interpreted from the symmetry 
argument and Pb-O bonding hybridization 37 . The first one is based 
on the fact that both the Pb0 12 and Ti0 6 polyhedra become increas- 
ingly distorted as temperature decreases. The distortion would intro- 
duce the lattice strain 21 , and as the temperature is raised and the 
polyhedra become more regular toward the Curie point, the unit-cell 
volume contracts continuously in response to the lattice strain 
relaxation. The second mechanism is complementary to the first 
one with more focus on the hybridization between the Pb 6s and O 
2p states. There are two distinct atomic positions for oxygen in the 
structure of tetragonal PbTi0 3 . Pb is strongly covalent with Oil, 
while the Pb-OI bond shows ionicity. The longer Pb-OII bonds 
alternate with shorter Pb-OII bonds along the c axis, as shown in 
Figure 4a. It is well known that as the unequal bond lengths become 
equal, their average value decreases, which, in the case of PbTi0 3 , 
would lead to the volume contraction as the tetragonal phase 
approaches the tetragonal-cubic phase transition temperature. 

To gain insights into chemical bonding and associated thermal 
behavior at high pressure we refined the atomic positions in tet- 
ragonal PbTi0 3 from the collected neutron diffraction data, as sum- 
marized in Table 1. Figure 4b shows the short and long Pb-OII bond 
lengths as functions of pressure and temperature. At all the tempera- 
tures the two Pb-OII bonds lengths are found to approach each other 
with increasing pressure, which is consistent with the prediction of 
density functional theory computations 38 . At any given pressure, 
Figure 4b reveals a similar tendency of bond equalization with 
increasing temperature. This is not unexpected because both heating 
at atmospheric pressure and compression at room temperature lead 
to the tetragonal- cubic phase transition. In addition, at all the tem- 
peratures, the average values of short and long Pb-OII bond lengths 
(as well as Ti-OI bonds) show similar trends of variation of cla ratios 
under pressure (Figure 3b). The observed trends for the bond-length 
equalization are also supported by high-pressure Raman spectro- 
scopic data 29 . It is known that the spontaneous polarization is pro- 
portional to the square root of (cla - 1). In Raman spectra this is 
reflected by the softening of the £(lTO) and A x (lTO) symmetry 
modes 29 with increasing pressure. By symmetry, the ions vibrate in 
the afr-basal plane in ^-symmetry modes, whereas they vibrate par- 
allel to the c-axis in A x symmetry modes. The softening of the £(lTO) 
and Ax ( 1TO) modes ultimately reflects the fact that pressure weakens 



Table 2 The mean thermal expansion coefficients of PbTiOs-based perovskite achieved by chemical 
at 0.3 GPa and 0.8 GPa are calculated from the interpolation of experimental data in Figure 3 

Compounds CTE(°C _1 ) 


and pressure-tuned routes. The values 
Temperature range 


PbTi0 3 , ambient pressure, Ref. 1 5 


-1.99X10" 5 


RT-763 K 


PbTi0 3 (0.3 GPa), this work 


0.47X 1 0" 5 


RT-600 K 


PbTi0 3 (0.8 GPa), this work 


0.31 xlO" 5 


400-700 K 


0.7PT-0.3BZT, Ref. 13 


-0.60x1 0" 5 


RT-773 K 


0.6PT-0.3BZT-0.1BF, Ref. 13 


-0.31 xlO" 5 


RT-773 K 


Fe[Co(CN) 6 ], Ref. 2 


-0.44x1 0" 5 


RT-773 K 


(l-x)PbTi0 3 -xBi(Ni 1 /2Ti 1/2 )0 3 , Ref. 14 


-0.88x1 0" 5 (x=0.1) 


RT-798 K(x=0.1) 




0.12X10" 5 (x=0.2) 


RT-798 K (x=0.2) 




0.85x1 0" 5 (x=0.3) 


RT-723 K (x=0.3) 
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Figure 3 | (a) A Evolution of lattice parameters a and c of tetragonal 
PbTi0 3 with increasing pressure at 300 K, 480 K, 580 K and 730 K. (3b), 
tetragonality (c/a) of PbTi0 3 as functions of temperature and pressure. 
The solid lines are drawn from linear fitting, and the determined pressure 
derivatives, d(c/a)/dP, are -0.013(2)/GPa, -0.0057(15)/GPa, 
-0.0039(16)/GPa and -0.00048(107)/GPa for 300 K, 480 K, 580 K and 
730 K, respectively. 

the bond strength and force constants responsible for the vibration 
frequency, which is directly correlated to the changes in ion positions 
and decrease in hybridization. Thus, based on the combined obser- 
vations of c/a ratios and equalization in Pb-OII bonds, the crossover 
behavior in volume isotherms is physically associated with tetragon- 
ality and hybridization between Pb and O atoms, which are both 
suppressed by pressure in tetragonal PbTi0 3 . 

Methods 

We conducted high-pressure/high-temperature neutron diffraction experiments on 
PbTi0 3 using a 500-ton, toroidal anvil press (TAP-98) at the beamline of High- 
Pressure- Preferred- Orientation (HIPPO), Los Alamos Neutron Science Center 
(LANSCE). Time-of- flight neutron data were collected using detector banks at a fixed 
angle of 26= ±90°. A PbTi0 3 powder sample was purchased from Alfa Aesar (99.9% 
metal based purity) and is phase-pure based on XRD. Two layers of PbTi0 3 and NaCl 
powders was first compressed into a cylindrical pallet of 5.4 mm in diameter and 
7 mm in length and was then loaded into a high-P- T ceramic cell assembly specially 
designed for TAP-98 (see Ref [39] for detail). NaCl served as the internal pressure 
standard, and the Decker's equation of state (EOS) was used to calculate the pressure. 
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Figure 4 | (a) A schematic illustration of crystal structure for tetragonal 
PbTi0 3 showing the alternative Pb-OII and Ti-OI bonds along the c axis. 
(4b), variation of short and long Pb-OII bond lengths as functions of 
temperature and pressure. Similar bond equalization tendency is also 
observed in Ti-OI bonds. The solid lines are drawn from linear fitting to 
guide eyes. 



Temperature was monitored by a W/5%Re-W26%Re thermocouple placed next to 
the sample and was stable within 25 K over a period of several hours of data 
acquisition time. To release the deviatoric stress built up during room-temperature 
compression on the polycrystalline sample, all high P- T neutron data were collected 
after the sample had been heated at 800 K for several minutes at each desired loading 
pressure. The high-temperature data at ambient pressure were collected with an ILL 
furnace; the sample was loaded into a vanadium can of 0.95 cm in diameter, which 
has low attenuation for neutrons and can hold temperatures up to 1500 K. The 
sample was heated under vacuum. The 144°-bank data obtained in ILL furnace and 
the 90° -bank data obtained in TAP-98 were used for determination of structural 
parameters of PbTi0 3 using the Rietveld method with the General Structure Analysis 
System (GSAS) 40 . 
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